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We use the bilayer twisted Hubbard model to study the low energy electronic structure of Sr3lr207. 
Sr3lr207 is suggested as a Mott insulator and described by the bilayer twisted pseudospin-1/2 model. 
The bilayer twists in Sr3lr207 bring the intrinsic anisotropy and the system is no longer SU(2) 
invariant. The anisotropy selects an easy c-axis coUinear antiferromagnetic magnetic structure at 
low temperatures which breaks the discrete Ising symmetry. The transverse spin excitations are 
gapped and the longitudinal one is weakly dispersive. We implement the bond-operator mean field 
calculation to describe the magnetic properties of the bilayer twisted pseudospin-1/2 model in good 
agreement with experimental measurements in Sr3lr2 07. The bilayer twisted Hubbard model is 
believed to be the good zeroth approximation for the low energy electronic structure of Sr3lr207. 



Introduction. The strong spin-orbit coupling in Ir^"*" 
(~ 0.4 eV) yields a narrow half-filled Jeff = 1/2 single 
band near the Fermi level in Ir-oxides and enhances the 
correlation effects of the one-site Coulomb repulsion (~ 2 
eV) [IHl] . It brings strongly correlated physics and makes 
Ir-oxides behave like Mott insulators, 1> iSfrlOj. The per- 
ovskite iridates are isostructural to the cuprate parent 
compounds and have been drawing increasing attention 
in recent years. Compared with the cuprates, the iridates 
are in the vicinity of the Mott transition^ [3 [TTl [T^ . 
With the help of the rapid-developing resonant x-ray 
magnetic scattering techniaue|13j. the low energy proper- 
ties of the iridates are well studied in the experiments [J 
[21 [51 ini [TIHTS] . The iridates provide the platform for in- 
vestigation of the low energy electronic structures of the 
Mott physics. The similarity and distinction between 
the cuprates and perovskite iridates can improve our un- 
derstanding of the Mott physics on the two dimensional 
square lattice. 

The perovskite iridates have Ir02 planes which deter- 
mine the low energy electronic properties, similar to the 
parent compound of cuprates. Ba2lr04 and Sr2lr04 have 
the single layer of Ir02 plane. As shown in Fig. [l] (a) 
and (b), the Ir-O-Ir bonds in the Ir02 plane are straight 
for Ba2lr04 and alternating twisted for Sr2l04, respec- 
tively. The low energy magnetic properties of Ba2lr04 
are the same as those of La2Cu04[71IH|- Sr2lr04 has the 
in-plane canted antiferromagnetic structure and shows 
ferromagnetism with large ferromagnetic moment due to 
the twists of Ir-O-Ir bonds [20]. The twists in Sr2lr04 
can be gauged away and the magnetic excitations are 
similar to Lao CuOa [20^22] . The activation energy gaps 
in electrical resistivity are estimated to be ~ 0.07 eV for 
Ba2lr04[7| and - 0.06 eV for Sr2lr04ll3], respectively. 
The perovskite Ir-oxides are in the vicinity of the Mott 
transition and the pressure can drive the systems cross 
the transition [3 [IH [H] . 

Sr3lr207 is the bilayer perovskite iridate with strongly 
coupled bilayers of Ir02 planes, stacked along the c- 
axis with pairs of adjacent bilayers offset by half a unit 
cell|24j. As shown in Fig. [l] (c), the Ir-O-Ir bonds in 



the bilayer Ir02 planes are twisted. The twisted direc- 
tions are alternating within each layer and also oppo- 
site between the two layers [IHl [23]. The optical gap in 
Sr3lr207 is finite although small[S], indicating the Mot- 
tness. Recently several measurements have been per- 
formed to explore the low temperature magnetic prop- 
erties of Br'^Ir^Ov fTlMTBl [^ . At low temperatures, the 
system has the long-range antiferromagnetic (AF) order 
as probed in the x-ray magnetic Bragg scattering piHTS] 
and the neutron scattering [5S]. The magnetic order has 
the easy c-axis collinear AF structure [HI [18], in the con- 
trast to the single layer Sr2lr04 with in-plane canted an- 
tiferromagnetic structure [rnj. The transverse spin wave 
excitations are gapped in Sr3lr207[Tn], in the contrast 
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FIG. 1. Schematic crystal structures of Ir02 planes in (a) 
Ba2lr04, (b) Sr2lr04 and (c) Sr3lr207. 
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to the gap less Sr2lr04[52]. There is a weak and broad 
feature above the transverse band in the inelastic x-ray 
scattering measurement in Sr-ilr907 |15l . The magnetic 
moment in Sr3lr207 is weaker than that in Sr9lrOA|18j. 
Increasing the temperature, the Bragg peaks of the AF 
order disappear at the magnetic transition temperature, 
Tat - 280 K[T3[25]. 

In this paper, we will use the bilayer twisted Hubbard 
model to describe the low energy properties of Sr3lr207. 
Due to the weak antiferromagnetism in Sr3lr207, the 
Mottness is questioned and the itinerant band magnet 
is suggested [m However, we insist on the Mott 

scenario for Sr3lr207. Although small, the optical gap 
in Sr3lr207 is definitely finite^. The large spin gap is 
not compatible with the weak magnetic moment in the 
itinerant band magnet [15] . Even if Sr3lr207 is the spin- 
density-wave band magnet, the magnetism has the same 
symmetry as that raised from the local pseudospins and 
the Mott scenario is still enlightening for the understand- 
ing of the magnetic structures. Due to the Mottness, 
we study the bilayer twisted pseudospin-1/2 model. The 
bilayer twists on Ir-O-Ir bonds in Sr3lr207 can not be 
gauged away and bring the intrinsic anisotropy. The 
anisotropy selects the easy c-axis antiferromagnetic struc- 
ture at low temperatures in which the discrete Ising sym- 
metry is broken. We implement the bond-operator mean 
field calculation to quantitatively describe the magnetic 
properties. The transverse spin wave excitations are 
gapped and the longitudinal one is weakly dispersive. We 
interpret the weak and broad feature above the transverse 
band as the weakly dispersive longitudinal spin wave ex- 
citations rather than the two-magnon scattering [15,. The 
interlayer coupling increases the itinerant of the electron 
and reduces the AF moment in Sr3lr2 07. On the mean 
field level, we ignore the quantum fluctuations and the 
spin gap and the magnetic moment reduction are signifi- 
cantly underestimated. For all that, however, the bilayer 
twisted Hubbard model still gives a good qualitative de- 
scription of the magnetic structures of Sr3lr207. 

Bilayer twisted Hubbard model for Srzlr2 O7. The crys- 
tal field in IrOg octahedra splits the d-orbitals of Ir"*"*" into 
Cg levels and t2g levels. The active t2g levels in Ir''+ has 
an efi^ective angular momentum 1 = 1, \lz ==0) = l-'^i^), 
\U = ±1) = -^{i\XZ) ± \YZ)) and is split further 
into Jeff = 3/2 and J^a = 1/2 by the on-site spin-orbit 
coupling. The ratios of the out-of-plane Ir-0 and in- 
plane Ir-0 bond lengths are 1.07 for Ba2lr04, 1.04 for 
Sr2lr04 and 1.02 for Sr3lr207l71 HH IM]. The elongation 
of the IrOe along c-axis is negligible and Ir^+ has the 
ideal Jcs =1/2 states as 

|J,^ff = +l/2) = -^(|0,t)-^/2| + l,;)), 

|J,^« = -l/2) = i=(|0,;)-V2|-l,t)). (1) 
As emphasized in Ref. 20, \XY), \Y Z) and \XZ) are 



the t2g orbitals defined in the local cubic axis of the IrOg 
octahera. 

In the Ir02 planes, the hoppings of Jeff =1/2 electrons 
between Ir sites are mediated by the oxygen 2p orbitals. 
For Ba2lr04 as shown in Fig. [l] (a), the Ir-O-Ir bonds 
are straight and the hopping term for J^s =1/2 electron 
on the Ir-O-Ir bond Tij is given as 

a 

where Cia is the electron operator on the z-th site with 
CT =t / i for Jgfj = ±1/2. For the straight Ir-O-Ir bonds 
the local cubic axis are also the global ones and the hop- 
ping parameter t can be specified as the real number. 

For Sr2lr04 as shown in Fig. [l] (b), the Ir-O-Ir bonds 
are twisted and the local cubic axis are alternatively ro- 
tated respect to the global axis. The twisted angle is 
9 ~ 11°. Here we still interpret Cia as the electron op- 
erator based on the local cubic axis [27]. On the twisted 
Ir-O-Ir bonds, the hopping term is now[31 [501 [2Z| 

/i*r' - E(^'""'"'^l^^- + 1^-^-), (2) 

with ei = ±1 for the sublattice and Eo- = ±1 for a =t / i. 
The hopping parameter is the complex number, te^*'^*'^"^. 
The twists of the single layer Sr2lr04 can be removed 
by proper site-dependent spin rotations and the twisted 
Hubbard model can be mapped to the SU(2)-invariant 
pseudospin-1/2 model[20j [21] which has the magnetic 
spectrum mimic that in the cuprate[22]. 

For Sr3lr207 as shown in Fig. [T] (c), we can read off 
the bilayer twisted Hubbard model similar to the above 
construction 

^ = E + E + E 4^^^l-AiCl^i^ (3) 

where {ij) denotes the nearest neighbor bonds in l-th 
{I = 1,2) layer. The hopping on the twisted bonds for 
each layer is 

M,--iE(^""'""'4.c/,.+h.c.), 

with ei = ±1 for the layer. The interlayer Ir-O-Ir bonds 
are straight in Sr3lr207 and the hopping on the interlayer 
rungs along c-axis is 

hi = -t± '^{cl^^C2^a + h.c). 

cr 

For Sr3lr207, the twist angle is ~ 12°. The hopping 
parameters {t and t±) are around 0.2 eV and the on- 
site repulsive interaction U takes the value around 2 eV 
analogous to those in Sr2lr04[2]. 

Based on the experiments (finite optical gap and large 
spin gap), we use the Mott scenario and obtain the bilayer 
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twisted pseudospin-1/2 model for Sr3lr207 



H = J2^[cos(26l)S,i • S,^ +2&in\e)S'{,S'{j 

{^3)1 

-e,ei sm{2e){Su x Sj^-) • z] + V Sh • S2.. (4) 



The super-exchange constants are given as J = At^ /U 
and Jj_ = 4i\/U, J ^ J^^c^m meV for SralrsOy. We 
can implement the site-dependent rotations and interpret 
the pseudospin in the global cubic axis 



'SI 
ySf: 



cos{e) e/e,sin(0) 0\ 
-e;e,sin(0) cos(0) | ^f, | , (5) 

01/ 



to remove the twists in each layer. Then we obtain the 
twists on the rung links as follows 



-2sin\6)Sl,S^^, 



e,sin(20)(Si,; x S2.) • z]. (6) 

When Jj^ = 0, the model ^ describes two decoupled 
SU(2)-invariant pscudospin-1/2 models as that in the sin- 
gle layer Sr2lr04 which has the canted antiferromagnetic 
structure and shows ferromagnetism with large ferromag- 
netic moment [m |20] . If we take the J± as the pertur- 
bation, we may anticipate the in-plane canted antifer- 
romagnetism in the bilayer Sr3lr207. The presence of 
the interlayer coupling prevents the full elimination of 
the twists in the bilayer twisted pseudospin model and 
brings the intrinsic anisotropy in Sr3lr2 07 which selects 
the easy c-axis AF magnetic structure different from that 
of Sr2lr04. When Jj^ ^ 0, however, the pseudospin align- 
ment along c-axis gains more energy from the interlayer 
couplings than the in-plane pseudospin alignment. The 
c-axis is the easy axis for the magnetic moment align- 
ment. Here we give alternative explanation of the easy 
c-axis different from that in Ref. [13]. Due to intrin- 
sic anisotropy of the bilayer twists in Sr3lr2 07, the low 
temperature AF order breaks the discrete Ising symme- 
try (not SU(2)) and the transverse spin wave excitations 
are gapped. The interlayer coupling J_l also introduces 
the quantum fluctuations and reduces the magnetic mo- 
ment. We will use the bond-operator mean field method 
to quantitatively describe these properties. 

Bond- operator mean field method. It has been proved 
that the bond-operator representation is a very conve- 
nient method to treat the bilayer spin svstem[28l |29] . 
For the rungs along c-axis, we introduce the bond 
operators '^B] (singlet boson s and triplet boson tx,y,z), 



NO) 



) - I m, \ty) 

410)- ' 

(Itt), 



1 



(I U) - I it)), \t.) ^ 410) = ^(tt 

4lO) = 75(1 tt) + I m and \t,} ^ 
yjd ti) + I it)) to describe the four spin states 
ti) I it, and I ii)). On every rung, we have the 



To simplify the calculations, we implement the mean 
field calculation at zero temperature and vary the in- 
terlayer couplings to tune the fluctuations. The singlet 
boson s condenses in all phase and the condensation of 
the triplet bosons describes the magnetic order. We start 
from the disorder state (large J± limit) and then decrease 
the fluctuations to see the magnetic instability. In the 
disorder state, we implement the approximation for the 
ground state, \G) = Hi ^i|0), and obtain the mean field 
Hamiltonian as 

i^MF = V j±tit,^ + ^ V(4.t,. + 4 A + h.C.) 
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(ij) 



^(cos(2^^)(4,^,,-4,^], + h.c.), 



(7) 



with a — x,y, z and t — x,y. The mean field Hamilto- 
nian can be diagonalized as i/MF = Eka ^kaCka^ka with 
the spin wave dispersions Wka = ^A^^ — F^^. Here we 

u A T I Jcos(2e) -p ,/cos(2e) A 
have Akr = Jj_ H 2^7k, Tkr = 2^7k, Akz = 

>/j_ + |7k, and T^^ = |7k with 7k = 2(cos(fc2,)-|-cos(fcj^)). 
Due to the intrinsic anisotropy of the twists, t^ has the 
smaller gap than tx,y at AF wave vector in the disordered 
state when Jl = 5 J as shown in Fig. [2] (a). Decreasing 
the fluctuations (J_l), t^ becomes gapless when Jl — 4 J 
as shown in Fig. [2](b). 

With further decreasing J_l, t^ condenses at AF wave 
vector and the system develops the easy c-axis collinear 
AF magnetic order. The AF order state is described as 
\G) —Y{i |0) with the new basis states, s\ = cos(x)s| - 
ti sin(x)i 



and i]. 



U sin(x)sl + cos(x)4.[3Q]. Now t, 
is the longitudinal spin wave excitation while t^ and ty 
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constraint s^s + Eq 4^a — \ [a 
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FIG. 2. Mean field dispersions of the bilayer Heisenberg spin 
waves for different values of ,J±: (a) J_i = 5 J, (b) Jj_ = 4 J, 
(c) Jx = J, (d) Jx = 0. The red (blue) line is for t^^y (t^) 
bosons. 
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are the transverse ones. The mean field Hamiltonian is 
now 



MF 



2q sin(2x)(cos(2x) - ^)(tl + Uz) 



i 



h.c. 



J 



f(cos(26') cos(2x)<Lijr 



2 ^ 



(cos(26') + i sm(26l) sin(2x)tL4r)) + ^-c]. (8) 



with T — x,y. Since \G) is the ground state, the coeffi- 
cient of the longitudinal spin excitation vanishes to give 
a self-consistent description, sin(2x)(cos(2x) — -jj) = 0. 
When Jj_ < 4J, the system is in the magnetic ordered 
state and cos(2x) = The spin wave excitations are 



:« = VKc - with = (2j + ^; 



J± cos(2e) 



7k, 



Tkr = 



J\ sin(2e 



-7k, Ak^ = 4J + W77k and 



Tkz = 32j7k- When < Jj_ < 4J, the transverse spin 
wave excitations have the gap 



A, 



V(4J + Ji)JiSin(0). 



(9) 



has the easy c-axis collinear long-range AF order, (ii) 
The transverse spin wave excitations are gapped and the 
longitudinal one is weakly dispersive, (iii) The magnetic 
moment is reduced by the interlayer couplings. The bi- 
layer twisted Hubbard model gives a good description of 
low energy electronic properties in Sr2lr307. 
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